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We extend the analysis of two-hadron fragmentation functions to the subleading twist, discussing 
also the issue of color gauge invariance. Our results can be used anywhere two unpolarized hadrons 
are semi-inclusively produced in the same fragmentation region, also at moderate values of the hard 
scale Q. Here, we consider the example of polarized deep-inelastic production of two hadrons and 
we give a complete list of cross sections and spin asymmetries up to subleading twist. Among the 
results, we highlight the possibility of extracting the transversity distribution with longitudinally 
polarized targets and also the twist-3 distribution e(x), which is related to the pion-nucleon a term 
and to the strangeness content of the nucleon. 
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I. INTRODUCTION 

(N 

' The study of the distribution of hadrons produced in the fragmentation of a quark offers the opportunity to 
understand the mechanism of hadronization as well as to extract information about the partonic structure of hadrons; 
C ■ both issues are a manifestation of confinement in QCD, a yet unexplained phenomenon. So far, parametrizations are 
available only for the distribution of the longitudinal momentum of only one of the final-state hadrons, the familiar 
unpolarized fragmentation function D\ (z) [lj, |2( • Clearly, most of the complexity of the fragmentation process lies 
unexplored. 

When the transverse momentum of one of the outgoing hadrons is measured, a new fragmentation function can 
be introduced relating the transverse polarization of the parent quark to the distribution of the produced hadron 
in the transverse direction [!|. This so-called Collins function acts as an analyzing power and it is perhaps the 
II 1 simplest observable that reveals the role of the quark's spin in the hadronization process. It also acts as a filter 
Q-i! to measure the still unknown distribution of transverse spin of quarks (transversity, for a review see Ref. Q) and 
the tensor charge of the hadron thereof However, the price to pay is the complete knowledge of the transverse 
dynamics of the detected leading hadron inside the jet. This creates problems both experimentally, as it is evident, 
and also theoretically, because the introduced dependence upon an intrinsic (nonperturbative) transverse momentum 
complicates the treatment of color gauge invariance 0, |(| and evolution equations 0, E lU • 

When two final-state hadrons are measured, in principle the number of variables doubles. For instance, it is 
possible to measure the relative transverse momentum of the pair, as well as its center-of-mass transverse momentum. 
Therefore, even after integrating upon the center-of-mass transverse momentum, a transverse vector is still available 
to establish a relation with the transverse polarization of the fragmenting quark 0, IT2I fl3| . 

Already from this intuitive discussion it is evident that two-hadron fragmentation functions can be important in 
studying spin effects in hadronization. They are perhaps more challenging to measure, inasmuch as they require the 
simultaneous detection of two hadrons inside the same jet. On the other side, the integration upon the center-of-mass 
transverse momentum removes the above mentioned difficulty about the evolution equations and it avoids, at least 
at leading twist, the potential loss of universality implied by a correct treatment of color gauge invariance 5, q, as it 
will also be discussed in Sec. IIIII 

Another class of functions that deserves much attention is that of polarized fragmentation functions. In this case, 
the spin of the final-state hadron is measured and its relation with the hadronization dynamics can be investigated. 
However, in general the spin of a final-state hadron can be analyzed only through the decay into two or more 
hadronic byproducts. In this sense, polarized fragmentation functions can be thought of as specific examples of multi- 
hadron fragmentation functions. For instance, the polarization of a vector meson (e.g., p°) is reflected in the angular 
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distribution of its decay products (e.g., 7r + 7r~). As a consequence, spin-1 polarized fragmentation functions [l4lll5llT6| 
correspond to the relative p- wave part of two-hadron fragmentation functions 17] . At present, however, the formalism 
of two-hadron fragmentation functions cannot comprise parity-violating decays, such as the extremely important case 
ofthe Abaryon [3 H H III EJ . 

Two-hadron fragmentation functions were first introduced in Ref. , but with no quark polarization. Extension 
of the original functions to include polarization effects (usually known as interference fragmentation functions) were 
studied in Refs. fTTllTp.123 . |2f| . The complete leading-twist analysis has been carried out in Ref. 01 an d employed in 
semi- inclusive DIS \26L |2JJ and electron-positron annihilation [28( . Positivity bounds and the expansion in the partial 
wave of the two hadrons were presented in Ref. 0] . Very recently, a study of collinear fragmentation into two hadrons 
has been performed [29|. demonstrating the factorization of two-hadron fragmentation functions at next-to- leading 
order in as and calculating their evolution, originally studied in Ref. |3Cj . In this article, we are going to extend the 
existing treatment to the subleading-twist level, but integrating upon the transverse momentum. The way we proceed 
is very similar to what was done in Ref. [3l|. for one-hadron production (see also Ref. |23), and in Ref. for the 
issue of color gauge invariance of the quark-quark correlator, even though we will only present results integrated upon 
the transverse momentum. The extension to the subleading twist is an important step not only from a formal point of 
view, but also because the measurement of two-hadron leptoproduction can be attempted in experiments at moderate 
Q 2 , where subleading-twist contributions should not be neglected. 

The paper is organized as follows. In Sec. [HI we will briefly review the kinematics for the semi-inclusive production 
of two unpolarized hadrons inside the same current jet. In Sec. IIIII we present the complete twist analysis up to 
subleading order of the quark-quark and quark-gluon-quark fragmenting correlators, discussing also the issue of color 
gauge invariance and of partial- wave expansion. In Sec. II VI the explicit expression of the hadronic tensor for the 
semi-inclusive production of two unpolarized hadrons in Deep-Inelastic Scattering (DIS) is shown, including leading- 
and subleading-twist contributions. In Sec.[V] the corresponding cross sections and spin asymmetries are discussed 
for different polarization states of the beam and the target. Finally, in Sec. I VII some conclusions are drawn. 




FIG. 1: Quark-quark correlation function A for the fragmentation of a quark with momentum k into a pair of hadrons with 
total momentum Ph = Pi + Pi and relative momentum R = (Pi — P-2)j2. 



II. KINEMATICS 



The fragmentation process is schematically represented in Fig. ^ where a quark with mass m and momentum 
k fragments into two unpolarized hadrons with masses Mi, M2, and momenta P\, P%. We introduce the vectors 
Ph = P1+P2 and R = (Pi —P^)/^. Using two dimensionless light-like vectors n+ and n_ (satisfying n 2 + = n 2 _ = and 
n + -n^ = 1), we describe a 4- vector a as [a~,a + ,a T ] in terms of its light-cone components =a-n^ = (a°±a 3 )/\/2 
and a bidimcnsional vector a T . 

For later use, transverse projection operators can be defined as 



tiv uv {m u } 

9t — 9 n + n - 



(1) 



where the braces indicate symmetrization upon the included indices and e 0123 = 1. In general, the fragmentation 
process is described in the frame where the hadronic final state has no transverse component, i.e., the frame where 
P hT — 0. We will see in the following that in actual measurements the natural choice is different and the required 
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boost introduces effects that have to be consistently taken into account when extending the analysis to the subleading 
twist. 

We also define the variables z = P^ /k~, the light-cone fraction of fragmenting quark momentum carried by the 
hadron pair, and the variable £ = 2R~ / P f ~ , which describes how the total momentum of the pair is split into the two 
single hadrons. Therefore, the relevant momenta can be parametrized as 



Ph <k 2 + k 2 ) % 



2PZ 



2Pr 



(2) 



C n _ (Mi - M|) - §M£ 



2 p h - 



IP: 



, R 7 



where Mh is the pair invariant mass. Not all components of the 4-vectors are independent. In particular, we note 
that 



Rf = 
R 2 -* = 

P h -R = 
Ph-k = 
R-k = 



2 _M 2 + Mi Ml 



(1 °i 1 + C) A4 2 - (1 - QM? - (1 + C)M| 



Ml - Ml 



(3) 



+ z- 



k ~\~ k rp 



Mi 

2z 

(Ml-Ml)-tMl 
2z 



The positivity requirement Rj, > imposes the further constraint 



M 2 h > 



1 + C 



Ml 



i-C 



-Mi 



(4) 



Note that to avoid the introduction of a new hard scale in the process, all invariants listed above have to be small 
compared to the hard scale Q of the process (where Q 2 = —q 2 , with q the momentum transfer). 



III. THE FRAGMENTATION CORRELATOR UP TO SUBLEADING TWIST 

The soft processes underlying the fragmentation are symbolically represented by the shaded blob in Fig. ^ and are 
described in terms of hadronic matrix elements of nonlocal quark operators as 

A(k,P h ,R)= ^ J (^I e ifc -« (0\^)\P h ,R;X){X;P h ,R\^(0)\0) , (5) 



where ip is the quark field operator. The above correlator is not color gauge invariant, as the quark fields are evaluated 
at two light-front separated space-time points, and £. To restore color gauge invariance, the so-called gauge link 
operator must be included, 

U [0,£] = ^ ex P (^~ [9 J Q dw-Mw)^ i (6) 



where A is the gluon field with coupling constant <?, and V indicates a path-ordered exponential. It symbolically 
corresponds to attach all possible soft gluon lines to the soft blob of Fig. ^ and resum their contribution. As such, 
the corresponding diagrams will still be considered as tree-level contributions, since the coupling g can be reabsorbed 
in the definition of the correlator itself. 
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The quark line in the fragmentation correlator has to come from a hard process that determines a dominant light- 
like direction. For the final state in a semi-inclusive DIS process, the hard scale Q selects the n_ direction as the 
dominant one with respect to the transverse and n + ones, which are suppressed as O(Mh) and 0(1/Q), respectively. 
The integration upon the suppressed n + components of the momenta can be performed up to 0{\/Q) leading to 

A(z,kr,R) = ^ J ^xfe ik < (0\m\Ph,R;X)(X;P h ,R\m\0)\ i _ =0 . (7) 

Similarly, in the gauge link it was usually assumed that the A + component of the gluon field is suppressed and, by 
neglecting the A T component and by imposing the choice A~ = (the so-called light-cone gauge), the gauge link was 
reduced to unity. Recently, the problem of the evaluation of such operator and of the gauge-invariant description of 
the quark-quark correlator has been studied in Refs. We will address it following the analysis of Ref. Q 

for the case of semi-inclusive DIS; the results can be easily generalized to the case of electron-positron annihilation. 
The proof in Ref. relies on counting the powers in 1/Q of the product of the fragmenting quark propagator and 
of the various components of the gluon fields attached to the soft blob, retaining only the leading and subleading 
contributions; these arguments are independent of the hadronic final state, and they are valid also for the two-hadron 
production, as long as the additional vector R does not introduce any new hard scale. 



If 















FIG. 2: Link structure for the leading-twist color gauge invariant quark-quark correlator for the fragmentation of a quark into 
a pair of hadrons. 



It turns out that only a combined analysis of leading and subleading contributions involving both A~ and A T 
components of the gluon field leads to a color gauge invariant expression for the fragmentation correlator. Following 
Ref. |5| (and generalizing its notation to the case of two-hadron production), a color gauge invariant object is obtained 
at leading twist by connecting the and £ points along the + direction running through — oo and through the transverse 
directions at £ + = — oo (see Fig.[2]for a schematic picture of the link path), namely 



A^(z,k T ,R) 
where 



x 



(2tt) 3 



l^-oo.g m \Ph, Rl X) (X; P h , R\ m U* tOB] CT+^j |0) 



£-=0 



(8) 



( 

\ 
/ 

\ 



ig / dw + A (w) 



w T — b T — a T 



ig I dw T ■ A T (w) 



o + = b+= a+ 
u~ — b~ — a~ 



(9) 



Note that by reabsorbing the product of gauge links, U T U + , into a redefinition of the quark field ip, the quark- 
quark correlator of Eq. JHJ falls back into the expression of Eq. Q , but for the [— ] superscript specifying the gauge 
link direction. Therefore, it still leads to a semipositive definite matrix in Dirac space |8fi| and the probabilistic 
interpretation of its leading-twist projections can be retained. The dependence on the direction of the gauge link is 
due to the contribution of U T , i.e. of the transverse component of the gluon field at £ + = — oo, which plays a crucial 
role in T-odd effects since it introduces nontrivial phases in the scattering amplitude. The direction of the gauge link 
depends on the considered process, potentially posing a threat to the universality of the definition of the soft correlator. 
For example, when considering the e + -e~ annihilation into one pair of hadrons in the same jet the correlator of Eq. © 
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will depend on a gauge link running through £ + = +00, therefore displaying the [+] superscript. However, it has been 
explicitly shown that universality is preserved at the one- loop level [(|. Moreover, when integrating AH(z, h?, R) of 
Eq. ijHJ upon dkx, the displacement of the quark fields is confined to the light-cone -I- direction. Hence, the two gauge 
links U T and U + will merge into a single operator connecting the points and £ along a straight line: 



A(z,R) = z 2 J dk T A H (z,£ T ,i?) = z 2 J dk T A [+] {z,k T ,R) 



z 2 



t x I ^ eiK <°l l7 [o,e]^)l^.^^)(^;^.,-R|^(o)|o) 



(10) 



£-=i T =o 



Therefore, in principle fc T -integrated functions are insensitive to the link path and should represent universal functions. 
This is certainly true at leading twist, but it is still matter of debate at the subleading level because of the presence 
of fc T -weighted contributions [3(|, as we shall see below. Therefore, here in the following we will omit an explicit 
dependence on the gauge link path for transverse-momentum integrated quantities only when the issue is settled and 
commonly accepted. 

At subleading twist, both the combinations of the transverse components of the quark propagator with A~ , and 
of the n + projection of the quark propagator with A T , generate a color gauge invariant operator involving the field 
strength tensor . After the fc T -integration, this correlator reads 



A [ A ]a (z,R) = [ dzr—^—AUz,zi,R) 

J Z\ + 16 



dz 1 ^—z 2 J( — ^e ik <e ik ^-i) (11) 
^i+ie fx J 2vr 2tt 



U+ $] ^(£) U+ n] gG- a ( V ) \P h ,R; X)(X; P h , R\ ^(0) |0) 



where z\ — /k~^ . This correlator starts contributing at twist 3; it is considered a tree-level contribution, as already 
explained at the beginning of this Section. Again, by absorbing the gauge link U + in a redefinition of the quark 
field i/j, the well known expression of the quark-gluon-quark correlator is recovered |37l l38| . Introducing the covariant 
derivative i£> M (£) = + gA^{^) we can recover also the relation 



where 



A [ A ]a (z, R) = A%{z, R) - A [ 9 - ]a (z, R) , (12) 

A a D (z, R) = z 2 ^j^- e ik < (0| U+ x] m iD a (0 \P h , R; X) (X; P h , R\ ^(0) |0) \^ =q , (13) 
A l d ]a (z,R) = z 2 J dk T k c ^A^-\z,k T ,R) , (14) 

with the covariant derivative -D(£) acting on the left on the quark field i/>(£). Note that, after integrating upon k T , 
the term A^, becomes insensitive to the gauge link path. It is possible to relate the quark-gluon-quark correlator to 
the quark-quark one using the equation of motion of QCD, (i Ip— m)ijj = 0. Therefore, A^ ' Q (z, R) does not introduce 
any new fragmentation functions, but it turns out that it plays an essential role in ensuring electromagnetic gauge 
invariance up to subleading twist. 

In the following, both the quark-quark and quark-gluon-quark correlators of Eqs. JSJ) and (|11|) . respectively, will be 
analyzed in detail for the semi-inclusive two-hadron production, including the expansion in the partial waves of the 
pair. 
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The correlator A 



The most general parametrization of A^(k, Ph, R) in Eq. (JSJ, compatible with Hermiticity and parity invariancc, 
is given by 

AM (k, P h ,R) = M h Cp 1 1 + 4 ±] Tf h + Cf ] # + Cf ] 



.3 

r i±] r l±] r [± 

+ -jj^— ft + -jj- (J^H ft + -jj- <7^f h H 

[±] 



^-^^ r— v-,," (15) 



where the coefficients are real scalar functions of all the possible independent invariants, namely ft 2 , k ■ Ph, k ■ 
R, Ml, Ml, Ml. Integrating Eq. (JT5J upon the suppressed k + direction and, consequently, taking the light-like 
separation £~ = 0, we get for, e.g., the DIS process the following decomposition: 



AH( 2) £ T ,ij;) = _L J dk+A^(k,P h ,R) 



k-=p-/z 



F pa f! k 



h 



(16) 



16ttP- i M h M h 

+ H l - 4 + ]+H*^ [#r> Vt] 

+ 75 M h +G 75 M h 

The first group of terms inside braces represents the leading-twist contribution and includes the usual inter- 
ference fragmentation functions (IFF) discussed elsewhere [3, LO H3- They can be obtained by projecting out 
of Eq. (|16fl the usual Dirac structures T — 7 _ 7s, icr i_ 75, where i means a transverse component. The 
second group shows the 1/P^ ~ 1/Q-suppressed fragmentation functions that arise from the Dirac structures 
r = 1, 7 l , <j h , cr^ , 7*75, respectively. Note that the structures T = i7s, a l+ , give no contribution at this level. 
The functions H?' , H^, , H, H^, G^', G 1 -, are T-odd, while H?' , , E, H, H^, are chiral-odd. 

Because of the constraints imposed by kinematics and by the k + integration, the fragmentation functions in Eq. I|16|) 
actually depend on five variables, namely z, £, M 2 , k^, k T ■ R T , and they can generally be decomposed as 

Di(z, C, Ml, kl k T ■ R T ) = D\{z, C, Ml, kl (k T ■ R T f) + D°(z, £ M 2 h , (k T ■ R T f) , (17) 

M h 

and similarly for the other functions. Both D\, D°, are even functions of k T . 
By integrating Eq. lltj|) upon the transverse momentum k T , we get 

A(z,R) = z 2 J dk T A^(z,k T ,R) 

h I M h 2 



(18) 



lGnPh I M h 2 11 M h 

A l (z, C, Ml, cf> R ) + A 2 (z, C, Ml, 4>r) 



where 



D <} = D <i'e + D ±o{l) j ( 19 ) 
G <i _ G </ e + G _Lo(l) j 
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and each term now depends on z, £, Ml We define the moment of a fragmentation function as 

Ht (1) (z, C, Ml) = / *t ^2 ^ C, ^ 2 , fcr, fcr ■ 4) - (20) 

and similarly for the other fragmentation functions. The resulting functions H^, H, G* 5 , are still T-odd, while 
, E, H, are chiral-odd. For later convenience, the leading-twist contribution is indicated by Ai and the subleading- 
twist one by A 2 , respectively. 



B. The subleading- twist correlator 

As already anticipated above, the color gauge invariant correlator A^" of Eq. (fTTfl is suppressed by one power of 
1/Q with respect to the leading twist AH of Eq. JSJ. Therefore, it must be consistently included when extending 
the analysis to the subleading twist. For sake of simplicity, only the fc T -integrated result will be shown. Since 

the gauge links can be absorbed in a redefinition of the quark fields, both Ag and A [ g ]q in Eqs. (T3j! and |QJJ, 
respectively, can be worked out in a way similar to the one-hadron emission. By projecting out the usual Dirac 
structures T = 7~, 7^75, icr l_ 75, the following decomposition results, 

A [ A ]a (z, R) = A a D (z, C, Ml cj> R ) - A [ g ]a (z, C, Ml <f> R ) 

= 1&kz I X/h" ~MT 75 ^ (21) 

l a i- ■ rr<o(l) Rt $T / 



where the functions with tilde denote 



D< = - z D° (1) , 

1 M h 1 



(22) 



H = H + 2zH^ 



.(i) 



and are pure twist-3 fragmentation functions depending on z, £, Ml They all vanish in the Wandzura-Wilzcek 
approximation. 



C. Partial- wave expansion 



If the invariant mass Mh is not very large, the hadron pair can be assumed to be in a channel corresponding to a 
relative s- or p-wave. Consequently, two-hadron fragmentation functions can be decomposed in partial waves 171. In 
the center-of-mass (cm) frame of the two hadrons, the emission occurs back-to-back and the key variable is the angle 
9 between the directions of the emission and of Ph- The kinematics described in Sec. [H] can be easily adjusted to the 
cm frame of the two hadrons; the most important modifications are 

Rrp = Rsm9 , 

= 2M~ h - 2(M i + M ') + (M? Mi) 2 > (23) 

C= W h ( V M i - l^l 2 - V M l - l^l 2 - 2|i?l cos<?) , 

where the crucial remark is that £ is at most a linear polinomial in cos 9 with coefficients that depend only on invariant 
masses. This suggests that the dependence upon £ in the fragmentation functions should be conveniently replaced 
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by an expansion in the Legendre polynomials in cos 6* and, consequently, the cross section kept differential in dcosO. 
The Jacobian d(/d cos 9 — 2\R\/Mh can be absorbed in a redefinition of the fragmentation functions. 

The partial- wave expansion of the leading- twist fragmentation functions has been given in Ref. 17], namely 1 



Di -> f 1,00 + cos e + D liW i (3 cos 2 9 - 1) 
H?^H? ot + H< lt cos 9. 



(24) 



Extending the analysis to the subleading-twist functions is straightforward: 



H 



<3o(l) 



H^+H^ cose, 



l l,ot 

D* 3 — * 5^ + 5^ cos , 
G^-G^ +G? t cos9, 

H^H OQ + H ol cos9 + H u \ (3cos 2 - 1) , 
£^£ oo + £ oi cos0 + £ n i(3cos 2 0-l) . 



(25) 
(26) 
(27) 

(28) 
(29) 



p - 



— I ^ 








< . 







- k 




FIG. 3: Relevant diagrams at leading and subleading twist for the SIDIS of a lepton on a hadronic target with detection of 
two hadrons in the same current fragmentation region. The shaded blobs understand the contribution of all unsuppressed 
longitudinal gluons, while the gluon lines represent all possible contributions from transverse gluon fields (see text). 



At variance with Ref. Il7l . here we use lower-case indices for the polarization of the relative partial wave, in order to avoid confusion 
with the polarization state of the beam and/or the target in the expression of the cross section (see the following Sec. IVt. 
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IV. HADRONIC TENSOR FOR SEMI-INCLUSIVE LEPTOPRODUCTION 



When the semi-inclusive production of two hadrons happens via a DIS process, an electron with momentum I 
scatters off a target nucleon with mass M, polarization S and momentum P, via the exchange of a virtual photon 
with momentum transfer q = I — V . Inside the target, it is assumed that the photon hits a quark with momentum 
p, changing it to a state with momentum k = p + q before the fragmentation (see Fig. We define the variable 
x = p + 1 P + , which represents the light-cone fraction of the target momentum carried by the initial quark. As 
already anticipated in Sec.^ it is customary to consider the frame where all the hadronic systems have no transverse 



components, i.e. where P T = P, 



h T 



0, while the virtual photon has a nonvanishing component q T . A convenient 



parametrization for the momenta referred to the initial hadronic system is 



M 2 



2P+ 
p 2 + 



,P + ,0 



Pt 



2xP^ 



xP^ 



,Pt 



(30) 




FIG. 4: Kinematics for the SIDIS of the lepton I on a (un)polarized target leading to two hadrons inside the same current jet. 

However, when calculating the hadronic tensor (and, consequently, the cross section) it is more convenient to 
consider the frame where the z axis is antiparallel to the direction of the virtual photon momentum (see Fig0J . By 
denoting the momenta in this frame with the subscript j_, we have, therefore, P± = q±_ = and Ph± — —zq T . The 
difference between the t and the _l frames is a boost that introduces corrections suppressed as l/Q; therefore, it can 
be neglected at leading twist, but it must consistently be included when extending the analysis at the subleading 
twist. The boost amounts to the following modifications, 

- rij 1 ~y^-q£ = + ^ (p T - k T Y , 

< ~ n» , (31) 

a T ~ 9± a T v q~ a T ' It n + = a T± "I" ~q a T ' \Pt ~ K T )n + , 

where a T is a generic transverse 4- vector, are the light-like vector considered in the j_ frame, and the analogue of 
the transverse projection operators of Eq. are 



9T =9^ ~n { ^n'_ v} 



(32) 
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As an example of the difference between the t and the ± frames, in Fig. [SI we sketch the vectors Rt and Rt±- 
As expressed in Eq. IpHJl. the difference between the two vectors is of order l/Q (exagerated in the drawing). The 
difference between the angles 4> R and <fi RJ _ and between |.R T | and |-R TJ J is of order l/Q 2 , therefore it can be neglected 
in our analysis. 




FIG. 5: Description of the angles <f> R and 4> R± . 
The hadronic tensor, integrated upon the transverse cm momentum of the hadron pair, reads 



2MW tUJ = 32zTr 

- 32zTr 

- 32zTr 

- 32z Tr 

- 32z Tr 



dp T dk T (x, p T , S) 7^ A H (z, k T , R) 7 " 



r ^= la A(z, R) 7" 7° * l Z ]a \x, S) 7° 
Y ^= 7o Hx, S) r 7° A^ 1 " 1 {z, R) 7° 
la -^=rA [ X ]a (z,R)r^,S) 



(33) 



Each contribution corresponds to a specific class of diagrams in Fig. [3] For sake of simplicity, the blobs in the 
diagrams understand all the connected lines related to unsuppressed longitudinal gluons, namely the lower blob 
includes all lines with A + gluons and the upper blob all lines with A~ gluons. Therefore, the diagram in Fig. [3^ 
corresponds to the first term in Eq. (|33|l involving the leading-twist color gauge invariant correlators $ (which will be 
described in Sec. lIV AT) and A of Eq. I|10|) . It is important to perform the integration upon the transverse momenta only 
after including the effect of the boost in Eq. (|31|) . as it will turn out that the correlator A contains a l/Q-suppressed 
p T dependence. 

At subleading twist, the contribution of transverse gluons A T is symbolically indicated by a line attached to the 

lower blob, corresponding to the color gauge invariant quark-gluon-quark correlators (which will be described 

in Sec. lIV A*) l. or to the upper blob, corresponding to A^'" in Eq. (|llf) . Therefore, the second and the third terms in 
Eq. i.YMi correspond to diagrams in Figs. EJ 3 and|3J:, while the fourth and fifth ones to diagrams in Figs. and^, 
respectively. 

The correlators A and A^ have already been discussed in Sec. 1111 Al and ITTl Bl respectively. In the following, 
the missing terms will be described in detail, leading to the final expression of in terms of distribution and 

fragmentation functions. 
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A. The quark-quark correlators for the initial and the final states 

The color gauge invariant quark-quark correlator for the initial state, &(x, S), corresponding to the lower blob in 
Figs. H2i, and|5k, reads 

$(x,S) = J dp T $ [+] (x,p T ,S) 

= \ 4+ + Sl9x(x)j 5 4+ + hi(x)7n $±4+ 

+ j 6 ^ + 9t{x)15 $± + S L h L (x)j 5 4+ 4- 1 ( 34 ) 

+ e L( x hs ~ f T ( x ) £t 7 Q S±p + 'M x ) 4+ 4- 

= $ 1 (x,S)+$ 2 (x,S) , 

where Sl/± are the longitudinal/transverse components of the target polarization, respectively. The first group (^i) 
represents the contribution of the leading-twist distribution functions and it appears in the first, fourth, and fifth, 
terms of Eq. (|33|l . corresponding to the diagrams of Fig.|3ji,[3jl, andJS^, respectively. The other terms ($2) represent 
the contribution of the subleading-twist distribution functions, including also the p T -integrated T-odd functions 
h(x), f T (x), and e L (x), which are vanishing if the gauge link is the only source of the T-odd behaviour. The $2 
contributes only to the first term of Eq. 1(33}, corresponding to the diagram of Fig. Since the non-integrated 
(x,p T: S) involves scalar products of transverse vectors and commutators between transverse vectors and the 
light- like vector n + , it is easy to check that the boost transformations in Eq. (|31[1 do not add other subleading-twist 
terms and leave $(a;, S) unaltered. Moreover, because of the p T -integration the latter is insensitive to the direction 
of the link integration path. 

At subleading twist, also the quark-gluon-quark correlator ${^ a comes into play, appearing in the diagrams of 
Figs. 13b and |3J;. Similarly to the previous case, the redefinition of the quark fields including the gauge links and the 
integration upon p T allows to keep the same relation — — ^a"'" dictated by the QCD equations of motion 

as in the non color gauge invariant case. The quark-gluon-quark correlator can be parametrized as |5j 

<&W>, S) = M I (xg T {x) - ^ hx(x) - g$(xj) SI 75 4+ 

+ S L (xh L (x) - -rjgi{x) + 2h^ 1) (x) +ixe L (x)) \ j 5 j a r/i + 

,m \ ( 35 ) 

+ [jj fi(x) - xe(x) + 2ihf [1) (x) + ixh(x)j i r 4+ 



(flT 1 \ x )+ x fT( x )) e< T S ±fl 4-\ 



where, as usual, we define the moments 



7-2 



g${x) = / dp T ^^g 1T (x,p^,) , (36) 



and similarly for the other distribution functions. 

As for the fragmentation into two hadrons, the leading- twist correlator Ai of Eq. (|18fl occurs in the diagrams of 
Figs.OtjEbj and 13:, while the subleading-twist A2 occurs together with $1 only in the diagram of Fig.|2t- However, 
the boost transformation l|31|) induces two additional contributions to the correlator of Eq. l|16f) . that are suppressed as 
1/Q and, therefore, must be consistently included in the analysis at subleading twist before performing the integration 
upon k T . The final £ T -integrated result reads 



2 V ' ' 16nQ \ 1 /5 M h 1 M h 



(37) 



ft A 2 a (z, R) = ^ p% I D lla + W< ^L- [ # T , 7q ] + i^ H< \[4'_ ,4 + ] } ■ (38) 
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Such contributions appear in the first term of Eq. (|33|l . corresponding to the diagram of Fig. The former couples 
to $ij while A 2ct , due to the presence of pj,, couples to 



<z> [ + ]a ( Xl s)^ J dp TP %& + Hx,p T ,s) 



y { /# (x)efS Tf) 4 + + (x) SI 75 i + 



(39) 



-s L h { l ) L {x) Tol a 4 + -\h^> 1 



(!), 



To complete the picture about the fragmentation at subleading twist, the quark-gluon-quark correlator A^ of 
Eq. H21|) must be included in the fourth and fifth terms of Eq. I|33|) . corresponding to diagrams OJi and|U- 



B. The hadronic tensor 



Putting together in a consistent way all the contributions discussed above up to the subleading twist, we get for 
the hadronic tensor the following expression: 



2M 



16z 

47T 



-9Tfi Di + i S l 9i D x - HT± - e± - ±p tr ^ ^ 



2A/ h 



xQ 2t^e v l P R T±p ( M „ G< 

+ Sl Q VM h XhLH i +91 ~ 



Q 



2tW&/ D< + M 



i Sj 



2i^ei p R T ±n ( D 



Q 

2Mt^t v } p S 1 
Q 



M 



9i -tt xe T ) — i 

z M h 



x qt D\ H All — 

y M 1 z 



Q \M h 1 n z 



The leading-twist contribution in the above formula involves color-gauge invariant quantities that are independent 
from the properties of the gauge link; under the hypothesis of factorization, it represents a universal response. At 
subleading twist, the issue is still under debate (EEJUll- However, it is interesting to note that only fc T -integrated 
fragmentation functions appear via Ag~' a in Eq. I|21|). leading to the "tilde" functions of Eq. I|22|) . that might depend 
on the considered process. No p T -integrated distribution functions appear via the corresponding 3>a~ , because these 

contributions in Eq. i|3*5|) are exactly cancelled by the ones generated by coupling A 2a of Eq. (|3*S|) to ^g 1 ^" of Eq. ll3*9"|) 
in the first, second, and third contributions of Eq. (|33() (see also Eq. (64) of Rcf. Q). The net result is that the 
functions h(x), f T (x), and e L (x), are the only T-odd distributions in the hadronic tensor and, if not vanishing, they 
must be generated by a dynamical mechanism that has nothing to do with the sensitivity to the link path. 



CROSS SECTION AND SPIN ASYMMETRIES 



The cross section for SIDIS of polarized leptons off polarized hadronic targets with two unpolarized hadrons in the 
same current fragmentation region, reads 



d 7 cr 



g?£ dMi dcj) R dz dx dy dip s 



E 



9 ? 

& ye a 
32zQ 4 



L^2MWr , 



(41) 



where a is the fine structure constant, y = (E — E')/E is the fraction of beam energy transferred to the hadronic 
system and it is related to the lepton scattering angle in the target rest frame, (j) s is the azimuthal angle of the target 
polarization with respect to the scattering plane, <f) R is the azimuthal angle of the R T vector with respect to the 
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scattering plane, measured either around the Ph direction or around the z direction (see Fig. |SJ. The indicated sum 
runs over the quark and antiquark flavors a. The hadronic tensor W£ v of Eq. (|40|l is contracted with the lepton tensor 



O 2 
V 



-2A(y)g^ + AB{y)^t v + AB^i^x" + + V(y)t^x^ 

+2iAC(y)ef -i\W(y)^e?% , 
where A is the lepton helicity, x the spatial unit vector, v 1 = {rv\ + n '^)/V2, and 



A(y) = 



i i f 



B(y) = (1 - y) , 



C{y) = y (J - 1 



V(y) = 2 (2 - y) y/T^y , 



(42) 



(43) 



W(y) = 2y^/T^j. 

For convenience, in the following we will indicate the unpolarized or longitudinally polarized states of the beam with 
the labels O and L, respectively. Similarly, we will use the labels O, L, T, to indicate an unpolarized, longitudinally 
polarized, transversely polarized, target. We can then deduce the following list of cross sections 2 : 



d 7 cr, 



oo 



d 7 <J OL 

dJ(jQrp 



2irQ 2 y 



J2el{A(y)f 1 (x)D 1 (z,C,M l 



V(y) cos0j 



I A, 



- z fi(x) D<(z, C, Ml) + -^-x h(x) H?(z, C, M 2 h ) 



2irQ 2 y 



\Rn 



x h L {x) H? (z, C, Ml) + ~ 9l (x) G« (z, C, Ml) 



(44) 



(45) 



2nQ 2 y 



M h 



+ V(y) sm<p s — 



r2\ , \ R t\ 2 rr<o(l) 



h!(x) [-H(z, C, Ml) + LjZL. H? ° w (z, £ M%) )-— x f T (x) D 1 (z, C, M, 2 ) 



M 



Mh 



2 

d7(j Lo = ^Q^ x J2 e * w (y) sin ^ 



Q 



^- x e(x) H? (z, C, Ml) (x) G< (z, £ M 2 h ) 



(46) 
(47) 



d 7 <J LL = A S L e l\ C (V) 9i{x) D 1 (z, C, M 2 h ) 



- z g x (x) D< (z, C, Ml) - x e L (x) H? (z, C, M 2 ) 



• ^ x g T (x) D x (z, C, Ml) - - hi (x) E (z, C, Ml) 
Mh z 



(48) 



(49) 



In the above formula, we stress again that we have included also the contributions of the /c T -integrated T-odd 
distribution functions h(x),f T (x), and e L (x), which are vanishing if the gauge link is the only source of a T-odd 
behaviour. It would be interesting to experimentally check this feature. 



2 The distribution and fragmentation functions are understood to have a flavor index a. 
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Several useful spin asymmetries can also be built out of the previous formulae. In Eq. I|46(l for d 7 a OT , the transversity 
hi can be isolated at leading twist through the fragmentation function in a sin((f> R + cf> s ) spin asymmetry. This 
asymmetry has been already discussed in leading-order analyses 0, |23 an d seems very promising with respect to the 
Collins asymmetry, since it does not need to keep memory of the k T dependence but rather of the direction of R T . 

While data from purely transversely polarized targets are not yet available, the HERMES collaboration has per- 
formed spin asymmetry measurements with targets longitudinally polarized along the lepton beam [3^, 0, , hence 
with a polarization 3- vector S = (S x ,0, S z ) in the lepton scattering plane (<p s = 0) and with a transverse component 
S x with respect to the direction of the momentum transfer along z. Because of the kinematics setup, S x is suppressed 
by l/Q with respect to S z [42j. In the present case of detection of two hadrons in the same jet, therefore, both the 
leading- twist d 7 a OT and subleading- twist d 7 a OL of Eqs. I|46[) and (|45|l . respectively, should be consistently considered 
at the same time when looking for a sin0 fl asymmetry. However, d 7 a OT is considerably simpler than the correspond- 
ing cross section for one-hadron SIDIS, because the information about the transversity is not contaminated by other 
contributions, as it happens with the Collins and Sivers effects, respectively. Moreover, in the Wandzura-Wilzcek ap- 
proximation the fragmentation function vanishes inside d 7 a OL ; therefore, a sin cf> R spin asymmetry for two-hadron 
SIDIS in the HERMES kinematics would approximately lead to the product of the fragmentation function times 
the transversity h\ and the distribution h L , which is anyway related to h\ itself via a Wandzura-Wilzcek integral 
relation. 

Again, if we neglect G^, a sin^ spin asymmetry with polarized beam and unpolarized target would give access 
to the chiral-odd distribution e(x), always through the chiral-odd fragmentation function H± , as it is evident from 
inspection of d 7 a LO in Eq. (|T7|l . The function e(x) has recently attracted a lot of interest |43|], because it is directly 
related to the soft physics of chiral symmetry breaking |44|. Its first isoscalar Mellin moment gives the scalar form 
factor. Although this form factor (describing the elastic scattering off a spin-i target via the exchange of a spin-0 
particle) has not yet been measured, its value at t = —Q 2 = 0, the so-called a term, can be deduced by low-energy 
theorems from the experimental pion-nucleon scattering in the time-like region at the so-called Chen-Dashen point 
t = -Q 2 = 2ml, witn m TT the pi° n mass HjJ 0, 0] . Unexpectedly, the a term turns out very big (50-70 MeV) |4al49| 
with respect to the average value of available lattice calculations [5(| , suggesting that approximately 20% of the nucleon 
mass M could be due to the strange quark content of the nucleon. Therefore, having experimental access to e{x) is 
of great importance. This distribution could be extracted at subleading twist through the Collins function by a beam 
spin asymmetry in one-hadron SIDIS for longitudinally polarized beams and unpolarized targets 0, |^| , provided 
that the transverse momentum of the detected hadron is measured. This asymmetry contains another contribution 
that was neglected until recently |52j,|53j. Once again, the case of one-hadron SIDIS is complicated by the dependence 
upon the partonic transverse momentum. For the case of two-hadron SIDIS, it is possible to integrate upon the 
transverse total momentum of the pair and still build an azimuthal asymmetry using R T . In fact, Eq. I|47|) looks 
simpler than the corresponding one for the one-hadron case, and it could eventually represent the cleanest channel to 
look at in order to extract e(x). 

Finally, when expanding the fragmentation functions in partial waves and making the cross section differential in 
cos 9, the different dependence upon 9 allows to distinguish the contributions pertaining to pure s waves, pure p waves, 
and s-p interferences. For instance, by substituting Eqs. i|24[l and (|23[1 into Eq. (|46|) it is possible to check that the 
asymmetry will be dominated by an s-p interference fragmentation function at 9 = tt/2, and by a p-wave interference 
fragmentation function at 9 = tt/4. 

VI. CONCLUSIONS 

Fragmentation functions are universal, process-independent objects |54j containing a crucial information about the 
hadronization mechanism and, ultimately, about the confinement of partons inside hadrons. They appear in semi- 
inclusive processes such as, e.g., DIS or electron-positron annihilation, and they can act also as a sort of "analyzing 
power" for the pola rization state of the fragmenting quark 0, [Tsl Il8l 132^ . The typical example is the so-called 
Collins effect |3J, l24j relating the transverse polarization of the parent quark to the transverse-momentum dependent 
Collins function, that describes a (nonperturbative) azimuthal asymmetry in the distribution of the detected leading 
hadron. Two-hadron fra gme ntation functions can also be defined, among which the so-called interference fragmen- 
tation functions [TH Il2l Il3j| lead to interesting single-spin asymmetries even after integrating upon the transverse 
total momentum of the pair |27l |28| , thus avoiding the complications introduced by the intrinsic nonperturbative 
transverse-momentum dependence of the Collins function. 

In this paper, we have extended the analysis of two-hadron fragmentation functions to the subleading-twist level, 
discussing also the issue of color gauge invariance but eventually integrating upon the transverse total momentum 
of the pair. Our results are theoretically interesting because the absence of an intrinsic nonperturbative dependence 
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upon transverse momenta cancels, at leading twist, also any dependence upon the properties of the gauge link operator 
necessary to restore gauge invariance, allowing for a truly universal definition of these objects; a debate is still ongoing 
to check if this property holds true also at subleading twist 0, 0, . The extension to the subleading twist is also 
experimentally important, because it can represent a nonnegligible contribution when performing measurements at 
moderate Q 2 . 

We have analyzed both the quark-quark and the suppressed quark-gluon-quark correlator, relating the latter to the 
former by means ot the QCD equations of motion. We have presented the full decomposition up to the subleading- 
twist level of these correlators in terms of fragmentation functions integrated upon the intrinsic transverse momentum. 
As previously stressed, these functions are universal certainly at twist 2 and maybe also at twist 3. 

As an application of our results, we have calculated the hadronic tensor and the cross section for all possible combi- 
nations of polarization states of the beam-target system in the case of deep-inelastic semi-inclusive leptoproduction of 
two unpolarizcd hadrons, by integrating upon the two-hadron center-of-mass transverse momentum. Our results can 
be used to distinguish 1/Q-suppressed contributions in experimental measurements, in order to extract more clearly 
leading-twist contributions, or in order to study interesting subleading-twist terms. An example of the former case 
is the possibility of extracting the transversity distribution in spin asymmetries also with longitudinally polarized 
targets (as they have been measured at HERMES |4l| for the case of one-hadron production); an example of the 
latter is the possibility of extracting from beam-spin asymmetries (probably in the cleanest possible way [52, 03) the 
twist-3 chiral-odd distribution function e(x) [5l|. related to the mechanism of the spontaneous breaking of the QCD 
chiral symmetry and, ultimately, to the strange-quark content of the nucleon |44j . 

As a last step, we have performed a partial-wave expansion of leading- and sublcading-twist two-hadron fragmen- 
tation functions, in order to distinguish the interference coming from the s-s, p-p, and s-p channels in the relative 
partial wave of the hadron pairs. Each component carries information on different mechanisms, such as the polariza- 
tion transfer to spin-1 resonances (for p-p interference) or T-odd effects from different kinds of final-state interactions. 
Therefore, extracting this information from data would allow for the exploration of different aspects of the physics of 
the fragmentation process. 
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